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Wo determine upper limits on the dark matter (DM) self-annihilation cross section for scenarios in 
which annihilation leads to the production of electron positron pairs. In the Galactic centre (GC), 
relativistic electrons and positrons produce a radio flux via synchroton emission, and a gamma ray 
flux via bremsstrahlung and inverse Compton scattering. On the basis of archival, interferometric 
and single-dish radio data, we have determined the radio spectrum of an elliptical region, around the 
Galactic centre of extent 3° semi-major axis (along the Galactic plane) and 1° semi-minor axis and a 
second, rectangular region, also centered on the GC, of extent 1.6° x 0.6°. The radio spectra of both 
regions are non-thermal over the range of frequencies for which we have data: 74 MHz - 10 GHz. 
Wc also consider gamma-ray data covering the same region from the EGRET instrument (about 
GeV) and from HESS (around TeV). We show how the combination of these data can be used to 
place robust constraints on DM annihilation scenarios, in a way which is relatively insensitive to 
assumptions about the magnetic field amplitude in this region. Our results are approximately an 
order of magnitude more constraining than existing Galactic centre radio and gamma ray limits. 
For a DM mass of = 10 GeV, and an NEW profile, we find {(Tav) < few x 10"^®cm^s"^ 

PACS numbers: 95.35.-|-d, 95.85.Bh, 98.70.Vc 



I. INTRODUCTION 

It is a remarkable fact that the identity of most of the 
matter in the Universe is unknown. An abimdancc! of 
observational evidence allows us to infer the existence of 
dark matter [1-3] via its gravitational influence. However 
we have as yet no direct detection and very little infor- 
mation about its corpuscular properties. In this paper, 
we focus on one of the fundamental particle properties of 
dark matter: its self annihilation cross section. 

If the dark matter were in thermal equilibrium in the 
early Universe, the annihilation cross section would de- 
termine the relic density in the Universe today. A ve- 
locity averaged annihilation cross section of {crAv)th ~ 
few X 10~^® cm^s~^ is required in order to produce the 
relic abundance of ^^dm ~ 0.2. However, it is not nec- 
essary for DM to be a thermal relic, in which case both 
(significantly) larger or smaller cross sections are possi- 
ble. Regardless of the thermal history, (ctav) controls 
the annihilation rate in DM halos in the Universe today, 
thus determining the size of detectable signals emanat- 
ing from regions of DM concentration. It is these anni- 
hilation fluxes that potentially permit the technique of 
indirect detection of dark matter. 

Recent cosmic ray positron/electron data from a num- 
ber of experiments have led to much excitement about 
indirect DM detection. Anomalies in the cosmic ray 
positron spectra have been reported by the PAMELA, 



Fermi and HESS experiments, implying an apparent ex- 
cess of positrons beyond those due to conventional astro- 
physical processes. PAMELA [4] has observed a rise in 
the e+/ {e~ +e~^) flux at energies above approximately 10 
GeV, while recent data from Fermi LAT [5] and HESS [6] 
show an excess in the (e^ + e+) flux up to and beyond 
1 TeV, respectively. (ATIC [7] and PPB-BETS [8] ob- 
served a similar excess, however, with considerably higher 
uncertainty than Fermi.) 

The explanation of these positron excesses is far from 
clear, and, in light of these, researchers have been moti- 
vated to examine or re-examine conventional cosmic ray 
interaction [9] and propagation models [10, 11], accelera- 
tion of e+e" in cosmic ray sources [12-14], electron and 
positron emission from supernova remnants [15, 16], and 
the production of positrons by pulsars [17-24]. As an 
alternative to conventional astrophysical mechanisms, it 
has also been speculated that DM annihilation or decay 
in the Milky Way may be the source of the excess elec- 
trons and positrons. 

Many authors have proposed models in which elec- 
tron and positron fluxes arise from DM annihilation [25- 
43] or decay [44-53], either directly or indirectly. (See 
also, Ref. [54] and references therein.) Note that, in 
contrast to the positron data, the PAMELA antipro- 
ton/proton observations do not indicate an anomalous 
contribution [57].^ Therefore, DM models which fea- 
ture signiflcant hadronic annihilation modes are con- 
strained, while leptonic channels are preferred. How- 
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^ For completeness, wc mention that both the reality of an anomaly 
in the PAMELA e+/(e^ + e+) flux and the absence of such in 
the anti-proton flux have been questioned [55, 56]. 
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ever, in order to account for the observed positron spec- 
tra, {(tav) must be significantly enhanced above the ex- 
pected thermal relic value. Such an enhancement could 
be of astrophysical origin, e.g., a boost due to signifi- 
cant substructure throughout the halo or a local clump 
of dark matter [58, 59]. Alternatively, it could be of 
particle physics origin, e.g., a Breit-Wigner resonant en- 
hancement [60-63] or the Sommerfeld effect (see, e.g., 
Refs. [25, 26, 64, 65]) in which low velocity annihilation 
(i.e. in Galactic halos) is enhanced while early Universe 
freeze-out is affected to a lesser degree [66-69]. 

Several techniques may be used to constrain the pro- 
duction of 6+ and e~ within Galactic halos, all of 
which rely on other accompanying observational signals. 
Charged particle production in a halo is necessarily ac- 
companied by photon signals, including gamma rays, X- 
rays, microwaves and radio waves. These signals are pro- 
duced via the various energy loss processes that charged 
particles undergo, examples of which include synchrotron 
emission in Galactic magnetic fields, inverse Compton 
scattering of electrons from interstellar radiation field, 
and bremsstrahlung. Charged particle production is also 
accompanied by internal bremsstrahlung radiation [70 
74] , which is an electromagnetic radiative correction, and 
is not due to interaction in a medium. 

In this work, we shall use the electron energy loss 
processes to derive constraints on {aAv). Annihilation 
channels we consider include both direct production of 
monoenergetic , and channels in which a spectrum of 
secondary are generated via decays of primary annihi- 
lation products, such as qq. Note that, given their identi- 
cal distributions and radiative signatures (at the energies 
under consideration), we shall imply both electrons and 
positron when we refer to "electrons" in this work, unless 
explicitly noted otherwise. Our analysis is distinguished 
from previous work in this area by the use of a new syn- 
thesis [75] of existing Galactic center radio data which 
allows us to derive particularly sensitive constraints. In 
addition, we make a careful study of the interplay of var- 
ious energy loss processes, and show that a combination 
of different techniques leave our final bounds relatively 
insensitive to uncertainties in magnetic field amplitude. 

A number of astrophysical uncertainties enter our cal- 
culations. As with all indirect detection analyses, we 
have a sensitivity to the assumed dark matter halo den- 
sity profile. In addition we are subject to uncertainties 
in the Galactic magnetic field amplitude, the background 
light density, and the gas density. Variation in the as- 
sumed magnitudes of these parameters alters the propor- 
tion of electron cooling that takes place via the various 
energy loss processes. Higher magnetic field strengths, 
for example, lead to greater synchrotron losses, while 
higher background light density leads to more inverse 
Compton scattering. However, while the astrophysical 
assumptions control the relative importance of the vari- 
ous electron energy loss processes, the eventual total en- 
ergy loss is a constant. (Note that, as we show below, 
relativistic electrons in the GC lose their energy - via 



whatever exact combination of processes - before they 
are transported out of the region.) Therefore, while con- 
straints from a single process (say, synchrotron radiation) 
are individually quite uncertain, we show that the com- 
bined constraints derived using all energy loss processes 
are robust and feature only mild dependence on these 
astrophysical parameters. 

This paper is structured as follows: in section II we dis- 
cuss the radio and 7-ray data which provide the empirical 
constraints on the various DM models we test. We also 
describe the physical environment at the GC (magnetic 
field, ambient hydrogen number density) which deter- 
mines how the cooling and subsequent radiation from rel- 
ativistic electrons proceeds and the DM profiles we inves- 
tigate. Section HI describes the distribution of electrons 
putatively injected by DM annihilations and explains how 
this injection spectrum is shaped into a steady state dis- 
tribution by the various cooling processes. We also set 
out here our calculation of the radio and 7-ray emission 
by these steady state electron distributions. In section 
IV we compare predicted emission from DM annihilation 
(parameterized in terms of DM particle mass and veloc- 
ity averaged cross section) against empirical data and 
delineate the regions of parameter space we can thereby 
exclude. We also compare our results with bounds exist- 
ing in the literature. Section V contains our concluding 
remarks. 



II. ASTROPHYSICAL INPUTS 
A. Radio Data 

LaRosa et al. [76] observed a discrete but diffuse, non- 
thermal radio source (DNS) covering a roughly elliptical 
region around the Galactic centre (GC) of extent 3° semi- 
major axis (along the Galactic plane) and 1° semi-minor 
axis between 74 and 330 MHz. This angular region corre- 
sponds, more-or-less, to the usual definition of the Galac- 
tic nuclear bulge [77] . Subsequently, the work of Crocker 
et al. [75] assembled archival, interferometric and single- 
dish radio data that demonstrates that the DNS radio 
structure is evident at frequencies up to 10 GHz and has 
a non-thermal radio spectrum over the full 74 MHz - 
10 GHz range of the data (see Fig. 1 for the spectrum 
and refer to Table I for references and other pertinent 
information on each radio datum). It should be noted 
though that, significantly for the purposes of Crocker et 
al., the spectrum of the DNS can not be described as 
a pure power law, but rather exhibits a spectral down- 
break at about a GHz. For a full account of the radio 
data processing and analysis supporting this conclusion 
the reader should refer to [75], but we repeat a number 
of particularly relevant points below. 

We also consider the constraints that arise from radio 
spectrum of the smaller region defined by |Z| < 0.8° and 
\b\ < 0.3°. We refer to this as the HESS region, as it cor- 
responds to a region for which the gamma ray intensity 
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is reported by the HESS collaboration. 



V (GHz) 


Telescope 


Beam 


Flux density 


Error 


Ref. 


.074* 


VLA 


2' 


16,200 Jy 


1,000 Jy 


[76] 


.330''^ 


Green Bank 


39 


18,000 Jy 


5% 


[76] 


1.408 


Effelsberg 


9.4' 


7,300 Jy 


10% 


[78] 


2.417 


Parkes 


10.4' 


4,900 Jy 


6% 


[79] 


2.695 


Effelsberg 


4.3' 


4,400 Jy 


10% 


[80] 


10.29 


Nobeyama 


2.9' 


1,400 Jy 


7% 


[81] 



TABLE I: Surveys and Observational data used to derive the 
spectrum for the 6° x 2° region centred on the GC. Notes: 
*At 74 MHz the large-angle Galactic plane synchrotron back- 
ground/foreground flux contribution is not measured (and 
hence not accounted for) due to the interferometric nature 
of the VLA observation. ^Total flux supplied by Dr Crystal 
Brogan (private communication). 
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FIG. 1: Radio spectra of the DNS (blue, upper) and HESS 
(red, lower) regions. 



Radio data at 74 MHz [76] were obtained by the Very 
Large Array (i.e., an interferometer), while data at all 
higher frequencies come from single dish instruments. 
(Because of the problem of free-free absorption along the 
line-of-sight to the GC [76] we do not employ the 74 MHz 
datum as a constraint in our fitting below. At higher 
frequencies, because of the c>c v~^-^ dependence of the 
free-free optical depth, absorption by this process is not 
a significant factor.) The radio flux at every frequency 
under consideration has or might have contributions from 
i) synchrotron emission by cosmic ray electrons along 
the line of sight, i.e., both in front of and behind the 
GC along the Galactic plane; ii) synchrotron emission 
by cosmic ray electrons gyrating in the GC itself; iii) 
discrete, astrophysical sources in the GC and along the 
line-of-sight (these include optically thin and optically 
thick thermal bremsstrahlung emitters and those with 
non-thermal spectra); iv) diffuse, large scale HH regions 
(i.e., regions of ionized atomic hydrogen that produce 
thermal bremsstrahlung); and v) (possibly and as inves- 
tigated in this paper) relativistic electrons and positrons 
injected following annihilation or decay of putative DM 
particles. 

As detailed in Ref. [75], the contribution of discrete 
sources (taken to be emission on angular scales < 1.2°) 
can be measured in the single dish data by Fourier anal- 
ysis. (At 330 MHz the contribution of discrete sources 
was measured directly with the VLA [76]). This discrete 
source contribution is at a < 20% level at all frequencies 
under consideration and has been removed from every 
radio datum we use in our analysis. Otherwise, we re- 
move no other astrophysical contribution and, in consid- 
ering possible limits on synchrotron emission by relativis- 
tic electrons and positrons injected by DM annihilation 
or decay, we take the conservative approach [82] of re- 
quiring that predicted synchrotron emission be no larger 
than 3(7 above any radio datum. 



B. Gamma- Ray Data 

In addition to synchrotron radiation, bremsstrahlung 
and (sub-dominantly) inverse Compton (IC) emission 
will inescapably be generated by any relativistic electron 
population given the presence of background gas (mostly 
molecular hydrogen) and the GC's background light field. 
Electrons with GeV scale energies produce both brems- 
strahlung photons of approximately 100 MeV and syn- 
chrotron radiation at GHz frequencies (given the range 
of magnetic field amplitude that is physically plausible 
for the GC; see below). This means that 7-ray observa- 
tions by the EGRET instrument (onboard the Compton 
Gamma-Ray Observatory) are particularly complemen- 
tary to the radio data introduced above, as both signals 
are produced by electrons of the same energy. Finally, 
additional 7-ray photons are also generated from neutral 
meson decay in any scenario involving quark-pair produc- 
tion in DM annihilation and both qq and €^ e~ scenarios 
generate internal bremsstrahlung 7-rays. 

1. EGRET Data 

We claim a conservative upper limit to the integral 
gamma ray intensity from the DNS region. This is ob- 
tained from Figs. 2(c) and 2(d) of Hunter et al. [84] 
which, when combined, show a pedestal in the super-300 
MeV intensity of ~ 3 x lO""*^ cm~^ s~^ sr~^, averag- 
ing over \b\ < 2° and within the longitude defined by 
[/[ < 30°. Given that the DNS region is much smaller in 
extent (in longitude) than the pedestal, and that there 
are no structures evident in these longitude-dependent 
intensity plots on the scale of the DNS, we estimate an 
upper-limit to any exotic contribution to the intensity 
in the DNS region at the level of twice the largest ac- 
tual excursion in the intensity above the pedestal value, 
viz. 1 X 10"'' cm~^ sr"-'. We are confident that a 
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more detailed analysis would produce a more stringent 
constraint. In any case, results pertaining to the region 
from the Fermi (GLAST) telescope are eagerly awaited. 



2. HESS Data 

Similarly, brcmsstrahlung and IC emission by the 
higher energy members of the electron population will 
generate ^TeV gamma rays (as will internal bremsstrah- 
lung and, where applicable, neutral meson decay). The 
differential intensity in this energy range may be com- 
pared with the (conservative) upper limit obtained from 
observations by HESS [83] . This instrument has detected 
diffuse emission over the approximate energy range 300 
GeV 10 TeV and over the smaller solid angle defined 
by 1^1 < 0.8° and |6| < 0.3°. The TeV intensity recorded 
is 1.4 X 10~^° cm~^ eV~^ s"-*^ sr~^ with only limited 
and dimmer diffuse emission detected outside this re- 
gion but within the DNS field. The HESS data analysis 
does, however, include a background subtraction from 
a nearby region, and this must be accounted for in our 
analysis[74, 85]. Because of this background subtraction, 
in the case of an NEW profile, a DM scenario may pre- 
dict an absolute flux of 7-rays at TeV energies and within 
the HESS field up to 2 x the observed ~TeV intensity. 
For the isothermal profile;, on the other hand, because 
of the very flat distribution and the consequently very 
small relative difference in 7-ray intensity predicted for 
the HESS region and the background region, no compet- 
itive constraint is proffered by the ~TeV data at all (and 
we, therefore, do not plot TeV 7-ray constraint curve for 
this profile). 



C. Ambient Environment 

1. Gas Density 

Over the DNS region, the volumetric-average hydrogen 
number density can be calculated to be about 13 cm~^ on 
the basis of the data presented in Ref. [77]^. Much of this 
density can be ascribed to molecular hydrogen found at 
very high number density (> lO'^'^ cm""^) but small filling 
factor (< 0.01) in the cores of the unusually-dense, Galac- 
tic center molecular cloud population [77, 86] . Quite pos- 
sibly electrons impinging from outside the dense molec- 
ular material may be excluded from it [87-89]. Exclud- 
ing this very high density phase {uh > 10^ cm~^) but 
including hydrogen in relatively low density molecular 



^ Wc assume the DNS volume to be an elliptical spheroid with 
circular cross-section in the Galactic plane (i.e., its greatest line- 
of-sight extent is equal to its width along the plane) . This implies 
a total volume of ~ 3.0 x lO'"'^ cm^. 



hydrogen, atomic hydrogen (HI) and 'warm' ionized hy- 
drogen (H+) one finds a minimum (path-integrated) nn 
for a DNS electron of 2.7 cm"^ and ~ 23 cm^'\ We 
adopt these values in our calculations henceforth. (The 
putative 'very hot', plasma phase of the DNS, that would 
have; a filling factor of ~ 85%, would have lower number 
density than this. However, relativistic electrons would 
not be trapped in this phase [90] and, in any case, the 
existence of this phases has recently been thrown into 
considerable doubt: see below.) 

2. Galactic Magnetic Field 

The amplitude of the GC magnetic field remains un- 
certain by about two orders of magnitude. 

• On the basis of the 74 and 330 MHz observations of 
the DNS (nuclear bulge) introduced above, LaRosa 
et al. [76] use the equipartition argument to suggest 
that the field on DNS size scales is only around 
6 /iG very similar in amplitude to that typical 
for the Galactic disk - climbing to only about 10 
liG over the inner l.°6 along the Galactic plane. 
(Coincidcntally, the region over which the HESS 
collaboration measures a diffuse ~ TeV 7-ray in- 
tensity [83].) These smaller field amplitudes have 
recently received support from rotation measures 
to external galaxies close to the GC [91]. 

• A rather stronger field over the DNS region is sug- 
gested by the work of Spergel and Blitz [92] who 
note that the putative "very hot" (8 keV) X-ray 
emitting plasma found throughout the GC [93] and 
the unusually dense and turbulent molecular gas 
might be in pressure equilibrium (implying near 
equipartition between these ISM phases). In anal- 
ogy with conditions local to the solar system [94] , it 
might then be expected that the magnetic field also 
contribute a roughly equal pressure requiring it to 
have about 100 /iG amplitude. The existence of the 
very hot plasma has recently been cast into extreme 
doubt, however, with the apparently diffuse X-ray 
emission being, it is contested [95], essentially as- 
cribable to unresolved, point X-ray sources. This 
would essentially nullify the Spergel and Blitz ar- 
gument. 

• Finally, (coincidentally or not) over the HESS field, 
radio observations of the 'non thermal filaments' 
(synchrotron-emitting radio structures, unique to 
the inner Galaxy, that run essentially perpendicu- 
lar to the Galactic plane) suggest that the ambient 
field might be as strong as 1 mG [96-98]. 

The situation regarding the GC magnetic field is entirely 
unclear, therefore. 

Recently this argument has been joined by Crocker 
et al. [75]. On the basis of a simultaneous analysis of 
the higher- frequency part of the radio spectrum of the 



5 



DNS radio structure and ^ GcV 7-ray data covering the 
same region (i.e., the very data described above), these 
authors have recently claimed a probable value for the 
magnetic field in this region of 100 /zG with a lower limit 
(at 95% confidence) at 50 fiG. However, we cannot self- 
consistently assume this result holds for the purposes of 
the current paper. The analysis in Crocker et al. as- 
sumes that the observed radio emission from the DNS 
region can be attributed to "conventional" astrophysi- 
cal electrons governed by a pure power law at injection, 
and that any spectral features detected in the observed 
spectrum arise in the cooling of the injection spectrum 
(which is initially featureless) . If we allow for the possi- 
bility that DM makes an additional contribution to the 
electron injection spectrum, the kinematics of DM an- 
nihilation allows for the possibility of spectral features 
in the injection spectrum of the in situ electron popu- 
lation, so the assumption of Crocker et al. that cooling 
alone introduces spectral features does not necessarily 
hold. (E.g., one can find a fit acceptable at better than 
2(7 confidence to the radio spectrum invoking synchrotron 
emission by a population of relativistic electrons injected 
following the annihilation of dark matter particles with 
a mass ^ 30 GcV and assuming the quark hadronization 
spectrum given below. We certainly do not claim this 
as a detection!) This implies we cannot self-consistently 
adopt the magnetic field strength inferred in Ref. [75] in 
this work. 

Given the uncertainties explained above, in this paper 
we investigate DM constraints across the whole range of 
magnetic field amplitudes supported by the existing lit- 
erature, which we take to be 10 - 100 /xG for the DNS 
region and 10 - 1000 /xG for the HESS diffuse flux region 
defined above. However, as we demonstrate below, ra- 
dio and 7-ray data offer complimentary constraints over 
this range: radio data provide the more severe restric- 
tion on {(Tj^v) for strong fields, while 7-ray data are more 
constraining for weak fields (because, in the latter case, 
there is less synchrotron suppression of the high-energy 
electrons required for IC up-scattering of ambient light 
into ~GeV and ~TeV energy ranges). 



D. Dark matter distribution 

Wc employ the usual Navarro-Prank-White (NFW) 
distribution [99] 

= ,2 > (1) 

- (1-h ^) 

where = 0.572 GeV cm~^ and rh = 14 kpc. The above 
implies a volumetric-average DM density of 3.6 x 10"'^" 
eV cm~'^ over the DNS region and 1.2 x 10^^ eV cm~^ 
over the HESS region. We employ these volume-average 
quantities in our calculations (given we assume a one- 
zone model). 



We also investigate limits in the case that the DM dis- 
tribution follows a truncated isothermal profile. Here we 
assume the parametrization presented in Ref. [82] follow- 
ing Ref. [100]: 

Pir) = ,2 ' (2) 

where ps = 1.16 GeV cm~^ and Vg = 5 kpc. 

III. DARK MATTER SIGNALS 

A. Dark matter annihilation and primary electron 
spectrum 

In typical models with weakly-interacting massive par- 
ticles (WIMPs), the DM masses fall in the GeV-TcV 
range, while the DM annihilation cross section is given 
approximately by the thermal relic value of {aAv)th ~ 
few X lO^^^cm^s"-'^. However, the DM self annihilation 
cross section can be much smaller if coannihilations play a 
significant role in determining the WIMP freezeout den- 
sity, as occurs in many supersymmetric DM scenarios. 
Conversely, the Sommerfeld effect can boost the self an- 
nihilation cross section relevant for galactic halos in the 
Universe today, effectively resulting in values much larger 
than the thermal relic expectation. Moreover, if dark 
matter was populated via some non-thermal mechanism, 
{crAv)th is not relevant to the presently observed average 
abundance. It is thus appropriate to take (cr^w) to be 
a free parameter, spanning a wide range of possible val- 
ues. Likewise, we shall also take the DM mass to be a 
free parameter. (Indeed, masses above 1 TeV have been 
considered in light of recent cosmic ray positron results 
and data from Fermi LAT.) We shall consider two DM 
annihilation processes: 1. XX ~^ QQ and 2. XX ~^ e+e~. 

i- XX m 

Annihilation to qq is taken as a representative proto- 
type of many DM models, as the resulting electron energy 
spectrum is a good approximation to that for many an- 
nihilation channels. For example, the annihilation spec- 
trum for a specific supersymmetric WIMP candidate may 
be calculated with the DarkSUSY package [101]. How- 
ever, the annihilation channels xX ^ qq, ZZ,W~^W~ 
all lead to essentially the same electron energy spec- 
trum (while the electron energy spectrum resulting from 
XX ~^ T~^T~ is significantly different). Note that for anni- 
hilation to ZZ, W^W~ , and qq, hadronization or decay 
of the primary annihilation products leads, via charged 
pions and muon intermediaries, to positrons and elec- 
trons with a broad spectrum of energies. 

For simplicity, we shall assume full annihilation into 
qq. For this channel, e~ and (e"*") will be produced by 
decaying muons (anti-muons) produced in charged pion 
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decay. The resulting e+ + e~ spectrum can be written, 
following Borriello et al. [102], in a simple polynomial 
form of the ratio Ee/m^c^: 



diVe 



(3) 



where J = {— |, — |, 0, |, 2, 3} and the coefficients aj are 
listed in Table II taken from [102]. 

This analytical expression is a reasonable approxima- 
tion for Eg > 1 GeV. At low energies the analytical ex- 
pression from [102] has an asymptotic behaviour cx E~^'^ 
while the kinematics of the charged meson and muon 
decay chains mean that the actual spectrum has a low- 
energy cut-off. Because of the large magnetic fields we 
sample, our calculated radio flux from synchrotron emis- 
sion is potentially sensitive to this effect at the lowest 
frequencies considered (cf. [102]). We therefore explic- 
itly introduce a cut-off by patching on a numerical de- 
scription of the low-energy spectrum we have previously 
obtained [87] in the context of a numerical treatment of 
secondary electron production from hadronic cosmic ray 
interactions (where identical kinematical considerations 
come into play) to Borriello et al.'s distribution. 

TABLE II: Uj values from [102] 



coefficient 


numerical 


1-3/2 


0.456 




-5.37 


ao 


10.9 


ffll/2 


-6.77 


a2 


0.969 


aa 


-0.185 



2- XX e e 

We shall also consider direct annihilation to monoen- 
ergetic electrons, for which the electron energy spectrum 
per annihilation is simply given by 

^jE,)^2Sim^~E,). (4) 

Note that the injection electron spectra differ from the 
steady state electron spectra due to the action of electron 
cooling, as discussed below. 



B. Steady-state electron distribution 

The injection spectrum of electrons produced in the 
Galaxy, at position r, as a result of Dark Matter annihi- 
lation is 

n /V _ 

(5) 



2 \ / dEe 



10'° 



^ 10' 



10" 10" 10'° 10" 10'^ 10" 

FIG. 2: Cooling time-scale due to various processes affecting 
electrons and positrons in the GC environment as a function of 
energy. Curves are - solid: ionization, dashed: bremsstrah- 
lung; dotted: synchrotron for (upper) 10 and (lower) 100 
fiG; and dot-dashed: inverse Compton. 



In general, the injected electrons loose energy interact- 
ing with the interstellar medium and move away from 
the production site. Both these broad processes lead to 
modifications of the shape of the injection spectrum. If, 
however, the transport timescale is much longer than 
the cooling timescale for electrons in the relevant en- 
ergy range one is in the so-called "thick-target" regime in 
which case the steady-state (cooled) electron distribution 
can be written as 



dEe 



(i?e,r) 



fl-'dE',Q{E',,r) 
-dE,{Ee)/dt 



(6) 



where dEg{Ee)/dt is the cooling rate, resulting from the 
sum of several energy loss processes that affect electrons 
(described immediately below). We show in the Ap- 
pendix that we are justified in adopting the thick target 
regime for the GC environment. 



C. Electron cooling processes 

In the GC environment four cooling processes are po- 
tentially important for charged leptons with energies 
such that they radiate at radio and 7-ray wavelengths 
(electrons and positrons suffer essentially identical en- 
ergy losses over the relevant energy scale and are treated 
identically). These processes are dominantly (from low 
to high energy) ionization, bremsstrahlung, and syn- 
chrotron and/or IC emission. We show the cooling 
timescales, as a function of electron (or positron) energy, 
in Fig. 2 for plausible GC environment parameters. In 
Fig. 3 we show examples of the shape of the injected elec- 
tron (and positron) spectrum together with the steady- 
state electron distribution after cooling. 
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Injection vs Cooled Spectrum 




10'^ 10' 10" 10' 10" 10" 

Es/eV 



FIG. 3: Shape of the injected electron + positron distribution 
(green; dashed) and steady state distributions (for various 
magnetic field amplitudes - mG (yellow; solid, upper); 100 
/iG (purple; solid, center); and 10 nG (blue; solid, lower). 
The curves have been normalized to 1 (on an arbitrary scale) 
at an energy of 10 GeV. Note the low-energy, kinematic cut- 
off in the injected distribution. The plot assumes = 100 
GeV. 



2. Bremsstrahlung Calculation 

We use standard formulae (see Ref. [104], p. 84 et seq.) 
to calculate the bremsstrahlung emissivity of relativistic 
electrons (and positrons) in the Galactic center environ- 
ment. 



3. Inverse Compton Calculation 

For inverse Compton cooling and emission calculations 
we assume the full interstellar radiation field determined 
for the inner 500 pc of the Galaxy in [105]. This has 
an energy density around 19 eV cm""^. We calculate the 
Compton scattering emissivity of relativistic electrons us- 
ing the full Klein-Nishina cross section set out in section 
2.3.3 of Ref. [106]. For the minimal nn value we assume, 
the integrated energy flux (100 MeV-f ) from IC emission 
is about 70% as large as the bremsstrahlung energy flux. 



4- Gamma-rays from neutral meson decay 



D. Radiative Processes 

Of the cooling processes listed above, all except ion- 
ization lead to potentially-observable signatures at radio 
and 7-ray wavelengths. The amount of cooling that takes 
places via either radio or gamma ray emission - and the 
constraints derived on the basis of that channel - depend 
on assumptions made about the astrophysical inputs. We 
calculate separate constraints from radio and gamma ray 
emission, in each case choosing astrophysical inputs con- 
servatively. In particular, for the hydrogen density, uh, 
we conservatively choose a value which leads to small- 
est observational signal for each process considered. The 
magnetic field strength is left as a free parameter to be 
varied, and we shall demonstrate the our final results are 
relatively insensitive to this parameter. In contrast, we 
fix the ambient hydrogen number density to be the ap- 
propriate minimal value set out in §11 C 1 (given the con- 
sideration that most of the ambient gas is 'locked-up' in 
molecular cloud cores of very high density but extremely 
low filling factor). We describe our treatment of each ra- 
diative process below. Representative 7-ray spectra from 
each of the relevant process (for a = 1 TeV DM mass 
and (cTAw) = 10~^^ cm^ s~^, with a ]3orriello et al. [102] 
spectrum) are shown in fig. 4. 



1. Synchrotron 

We use standard formulae [103] to calculate the radio 
emissivity of relativistic electrons (and positrons) in the 
Galactic center magnetic field. 



For DM annihilation into qq as investigated by Bor- 
riello et al., a 7-ray signal from the decay of neutral 
mesons produced in quark hadronization is inescapable. 
The contribution of this process is often dominant for rel- 
evant energies and environmental parameters (see fig. 4) . 
We employ a parameterization of the numerical results of 
Bergstrom et al.[107] for the 7-ray spectra produced by 
quark hadronization following annihilation of 500 GeV 
neutralinos (displayed in the their figure 12 and param- 
eterized in terms of the scaling variable x = Ej/m^). A 
representative distribution is that following uu produc- 
tion: 

^ = exp(-6.8x) (7) 
ax 

We also assume a branching ratio of 100% into the max- 
imally massive quark pair allowable given the DM mass 
(which means either bb or tt over the range of we 
explore) . 



5. Gamma-rays from internal bremsstrahlung 

DM annihilation to any charged particles is necessar- 
ily accompanied by internal bremsstrahlung (IB) radia- 
tion [70-74]. This is an electromagnetic radiative cor- 
rection, and not to be confused with the regular brems- 
strahlung considered in this paper. IB produces hard 
gamma rays up to the dark matter mass, with an approx- 
imately model-independent spectrum, largely free from 
both particle and astrophysics uncertainties. For the 
(sub-dominant) contribution of internal bremsstrahlung 
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FIG. 4: Sample 7-ray spectra of the DNS region for a = 1 
TeV DM mass and (aAv) = 10"^^ cm^ s"\ with a Bor- 
riello et al. [102] spectrum, a plausible ambient field of 100 
/iG, and an ambient hydrogen number density of ~ 3 cm~^. 
Curves are: (blue, solid) inverse Compton; (purple, dashed) 
bremsstrahlung; (yellow, dotted) neutral meson decay (fol- 
lowing hadronization) ; and (green, dot-dash) internal brems- 
strahlung. 



we employ the parameterization of Refs. [71, 74]: 



ttE^ 



In 



(8) 



where s — Arn^ and s' 



4:m^{m^ — E). (This equation 
has to modified for the case of qq production by replacing 
the electron mass with the effective quark mass, and a 
with a, where Qq is the electric charge of quark q [72] .) 



IV. DM ANNIHILATION CONSTRAINTS 
A. Discussion of results 

Figures 5, 6, 7 and 8, display the bounds we obtain for 
the annihilation cross section, as a function of DM mass. 
We show results corresponding to our four main scenar- 
ios, namely, the XX ~^ QQ ^^id XX ~^ e+e^ processes, 
for both NFW and isothermal DM halo profiles. For 
each scenario, we show how the constraints depend upon 
the assumed magnetic field strength within the plausible 
range 10 - lOO^iG. 

Neglecting the ^TeV 7-ray constraints for the moment, 
in each scenario investigated 7-ray data is the most con- 
straining at low magnetic field amplitude, and radio the 
best at high. Adjusting the magnetic field amplitude al- 
ters the proportion of electron cooling that occurs via 
synchrotron emission. Thus, for low (high) B amplitude 
the radio flux is less (more) constraining. The total en- 
ergy loss of the electron population is controlled by a 
combination of the available cooling mechanisms, such 



that reducing the amount of cooling taking place by syn- 
chrotron emission leads to relatively more gamma rays 
emission via bremsstrahlung and inverse Compton scat- 
tering. (Note that any individual electron must eventu- 
ally loose all its energy via one mechanism or another, 
and it is impossible for this energy loss to occur without 
the production of detectable signals.) Therefore while the 
radio constraints, considered alone, are quite sensitive to 
the assumed magnetic field strength, our final results, 
being the strongest of either the radio or gamma con- 
straints at a given B amplitude, are much less sensitive 
to the assumed value of B. 

For the NFW profile, we find that the ~TeV 7-ray 
constraint invariably becomes the best for larger values 



of 



(As remarked above, because of the background 



subtraction procedure, TeV 7-ray data does not proffer a 
competitive constraint for relatively flat profiles like the 
isothermal). The '^TeV 7-ray constraint curve is highly 
structured as a result of it being a comparison between 
predicted and observed differential 7-ray intensities at 
the nine separate energies the HESS collaboration has 
listed starting at ~300 GeV in concert with the three 
(for e^e^) of four (for qq) 7-ray production mechanisms 
and their different thresholding effects with respect to the 
parameter. A somewhat similar effect - though less 
pronounced is see in the radio constraints. In contrast, 
the EGRET 7-ray constraints (taken above a single pho- 
ton energy of 300 MeV in the integral intensity) define a 
fairly smooth constraint for each scenario whereas . 

In Fig. 9 we compare the best constraint obtained 
(from any modality except for ~TeV 7-rays) for plau- 
sible magnetic field values, namely 10, 30, 100, 1000 
/iG. We see that the best constraint on (av) varies by 
a factor of < 10 as S is varied from 10-100 fiG. (This 
is to be compared with a larger variation in the radio 
constraint alone of ^ 30.) Therefore, when considered 
in totality, the 7-ray and radio data conspire to gener- 
ate constraints on DM mass/annihilation cross section 
that are roughly constant across the reasonable range of 
astrophysical parameters (in particular, magnetic field 
amplitude). For larger magnetic field amplitude (i.e. 
B = 1000/iG, which is at the upper end of the astrophys- 
ically plausible range) synchrotron emission dominates 
the cooling, and somewhat stronger constraints are ob- 
tained. This occurs because the field has become strong 
enough that the GHz range synchrotron-radiation is sam- 
pling lower energy (and therefore more numerous) elec- 
trons. 

Comparing the figures for the NFW and isothermal 
profiles (e.g. Fig. 5 with Fig. 7, or Fig. 6 with Fig. 8) it 
is evident that the exclusion curves for the NFW profile 
are always much more constraining than for the corre- 
sponding isothermal scenario. This is due to the high 
central DM density of the strongly peaked NFW profile, 
in comparison with the low central DM density of the 
fiatter isothermal profile 

Likewise, the radio observation of the HESS region 
(small solid angle at the GC) provide more stringent con- 
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(a) XX^ qq; NFW profile; B = lO^tG. 



(b) XX qq; NFW profile; B = 30/iG. 



(c) XX -s- qq; NFW profile; B = lOOAtG. 



FIG. 5: Exclusion contours for the DM annihilation cross section assuming a XX ~^ annihilation channel (with a Borriello 
et al. [102] spectrum) and an NFW DM profile, obtained from DNS radio data (red, solid), HESS region radio data (blue, 
dotted), DNS region ~300 MeV 7- ray data (green, dashed), and IfESS region ~TeV 7-ray data (orange, dot-dashed). Relevant 
7-ray production processes are neutral meson decay, bremsstrahlung, inverse Compton emission, and internal bremsstrahlung. 
Results are displayed for three Galactic magnetic field amplitudes within a plausible range (a) 10 /iG, (b) 30 ^G and (c) 100 




FIG. 6: Exclusion contours for the DM annihilation cross section, assuming a xX ~^ e^e~ annihilation channel (mono-energetic 
spectrum) and an NFW DM profile, obtained from DNS radio data (red, solid), HESS region radio data (blue, dotted), 
DNS region ~300 MeV 7-ray data (green, dashed), and HESS region ~TeV 7-ray data (orange, dot-dashed). Relevant 7-ray 
production processes are bremsstrahlung, inverse Compton emission, and internal bremsstrahlung. Results are displayed for 
Galactic magnetic field amplitudes (a) 10 ^G, (b) 30 ^.G and (c) 100 /iG. 



straints than the radio observations of the DNS region 
(larger solid angle) for the peaked NFW profile. For the 
isothermal profile, however, we have something of the op- 
posite situation: the radio bounds for the DNS regions 
are often more constraining than those for the HESS re- 
gion. This can be understood since, for the flat isother- 
mal profile, the average DM density is approximately the 
same for both regions, while the radio background for the 
DNS region is lower than for the HESS region. 

With the exception of the HESS '-^TeV 7-ray con- 
straint, the qq injection scenario tends to generate some- 
what tighter constraints than the monoenergetic e^e~ 



spectrum, for all values of larger than ^ 100 GeV. In 
the qq case, electrons and positrons are injected as secon- 
daries from meson decay following quark hadronization. 
Therefore, a single xX ~^ interaction produces more 
electrons and positrons (though of much lower energy) 
than for the monoenergetic xX ~^ e^e^ case, and thus 
generally tighter limits. In addition, the production of 
7-rays from neutral meson decay, leads to strong MeV 
gamma ray limits in the case of qq. At the lowest ener- 
gies (below ~ 100 GeV), the e^e~ radio limits surpass 
those for qq, due to the higher average energy of the in- 
jection spectrum. 
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i'^) XX ~^ 11\ Isothermal profile; B = lO^G. (b) XX ~^ 1Q\ Isothermal profile; B = 30/iG. 



(c) XX ~^ Q<lt Isothermal profile; 
B = lOOAtG. 



FIG. 7: Exclusion contours for the DM annihilation cross section assuming a xX ~^ W annihilation channel (with a Borriello 
at al. [102] spectrum) and an Isothermal DM profile, obtained from DNS radio data (red, solid), HESS region radio data 
(blue, dotted), and ~ 300 MeV 7-ray data (green, dashed). [For such a flat profile the ~TeV data do not offer a competitive 
constraint: see the text.] Results are displayed for Galactic magnetic field amplitudes (a) 10 fiG, (b) 30 fiG and (c) 100 /iG. 




FIG. 8: Exclusion contours for the DM annihilation cross section assuming a xX ~^ e^e~ annihilation channel (mono-energetic 
spectrum) and an Isothermal DM profile, obtained from DNS radio data (red, solid), HESS region radio data (blue, dotted), 
and ~ 300 MeV 7-ray data (green, dashed). [For such a flat profile the ~TeV data do not offer a competitive constraint: see 
the text.] Results are displayed for Galactic magnetic field amplitudes (a) 10 fj,G, (b) 30 /xG and (c) 100 /xG. 



B. Comparison with previous results 

A number of previous works have considered DM anni- 
hilation corss section constraints obtained from Galactic 
radio and 7-ray data. While our analysis is complemen- 
tary to these prior studies, distinguished by the use of dif- 
ferent GC radio data, it is useful to compare the strength 
of the limits obtained on (av) : 

1. Of recent work, the most immediately comparable 
to ours - as discussed above - is that of Borriello et 
al. [102] whose NFW profile and electron injection 
spectrum we adopt (for the xX ~^ process). Our 
results are at least one order of magnitude and up 
to two orders of magnitude better than those ob- 



tained by these authors. The constraints of Bor- 
riello et al. are obtained from consideration of all- 
sky radio data at frequencies between 100 MHz and 
23 GHz. These authors, however, specifically ex- 
clude a 15° X 15° patch of sky centered on the 
GC (because their calculations are performed in the 
limit that only synchrotron cooling is important in 
shaping the steady state electron distribution and 
- as discussed at length above - in the higher gas 
densities near the GC, this will no longer be true 
given the consequent importance of ionization and 
bremsstrahlung) . 

2. Hooper [108] has considered constraints obtainable 
from WMAP data [109] at 22, 33, and 61 GHz. 
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FIG. 9; Plot of the best constraint we find (from any modality 
except fot ~TeV 7-rays) over plausible magnetic field values, 
assuming a XX ~^ annihilation channel and an NFW DM 
profile. Curves are; green - 7-ray constraints for (dashed) 30 
fiG and (solid) 10 ^G; and blue - radio constraints from the 
HESS region for (solid) 100 fiG field and (dashed) 1 niG field. 
Also shown for comparison is the best constraint obtained by 
Borriello et al. [102] for this scenario, purple, dotted. 



His analysis relies on astrophysical foreground- 
subtracted data that putatively reveals a residue, 
unexplained "haze" [110] at all three WMAP fre- 
quencies on scales out to about 20° from the GC 
(and absent elsewhere). The Hooper analysis, then, 
is intrinsically less conservative than the radio anal- 
ysis presented here where we do not attempt such 
subtraction (except for discrete sources in the field 
which can be unambiguously identified on angular 
size grounds). While Hooper's constraints are par- 
ticularly strong, some of the astrophysical assump- 
tions adopted may be too optimistic. In particular 
he assumes a constant 10 fiG field and constant 
5 eV cm~^ interstellar radiation field energy den- 
sity over the 30° angular region he investigates, the 
former being rather stronger and the latter rather 
weaker (at least towards the GC) than assumed 
by other authors [102, 111]. Grajek et al. [Ill] 
consider the same (foreground-subtracted) WMAP 
data investigated by Hooper. Assuming a rather 
weaker (and spatially-varying) magnetic field and 
a rather smaller ratio Ub/{Ub + C^rad) ^ 0.1, these 
authors find correspondingly weaker constraints on 
the annihilation cross section, which we surpass. 
(Note that the W^W~ annihilation channel con- 
sidered by Hooper and Grajek results in a very sim- 
ilar electron injection spectrum to the qq channel 



we consider.) 

3. Finally, Bertone et a/. [82] and Regis and Ullio [112] 
both also consider constraints arising from multi- 
wavelength observations of the GC, including ra- 
dio and 7-rays. These papers tend to focus on ob- 
servations on much smaller angular scales around 
the super massive black at the Galactic dynamical 
center than investigated here. While some of the 
constraints (and projected constraints) derived are 
quite strong, on these small scales it is necessary 
to assume a particular model for the evolution of 
the magnetic field intensity (and the matter density 
controlling bremsstrahlung and ionization cooling) . 
Such constraints are also particularly sensitive to 
the DM density at extremely small radii, which is 
highly uncertain. 

We now compare with constraints derived using other 
techniques. For hadronic annihilation modes, constraints 
arising from the antiproton observations have been shown 
to be of comparable sensitivity to existing Galactic centre 
gamma ray and radio limits, e.g. Refs. [113, 114], for 
some models. For the monoenergetic e"'"e~ annihilation 
channel, pure IB constraints were derived in Rcf. [74]. 
For the NFW profile, the IB constraints are very strong 
for 10'^ < < 10"' (these arise from HESS gamma ray 
data near the GC) whereas our new radio/gamma ray 
bounds are better at lower mass. (For an isothermal 
profile, the IB constraints arising from the HESS data 
would be considerable weaker.) Note that for masses 
above ~ 10 TeV, bounds on the total annihilation cross 
section to all final states, derived using neutrino data, 
become most constraining [115-117]. For higher masses, 
™x > 10^, TeV the unitarity bound on the total cross 
section becomes the most restrictive constraint [118]. 

Turning now to non-Galactic based techniques, a num- 
ber of authors have recently considered constraints aris- 
ing from annihilation at high redshift. For example, 
Refs. [119-123] consider the effect of energy deposited 
by DM annihilation during the the reionization and re- 
combination epochs. (Note that these techniques are 
complementary to those we consider here, being subject 
to quite different assumptions and systematics.) For an 
NFW profile and an e+e" annihilation channel the reion- 
ization limits are broadly comparable to the constraints 
we have derived here, with the boost factor required to 
be < 10 at - 100 GeV. 



V. SUMMARY AND CONCLUSIONS 

We have considered DM annihilation in the Galactic 
center, for scenarios in which annihilation leads to the 
production of relativistic electrons and positrons. Us- 
ing signals arising from electron energy loss processes, 
namely synchrotron emission, bremsstrahlung and in- 
verse Compton scattering (and, where appropriate, those 
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due to neutral meson decay and internal bremsstrah- 
lung) , we have derived robust constraints on the velocity- 
averaged DM self annihilation cross section {cv) . The 
processes considered were XX e+e~ which produces 
monoenergetic e^, and XX ~^ QQi in which a spectrum of 
are generated via the charged pion decay chain. (The 
constraints for the later channel are expected to be very 
similar to those for other interesting final states, such as 
W+W- and ZZ.) 

We have demonstrated that a combination of radio 
and gamma ray bounds is relatively insensitive to the as- 
sumed magnetic field amplitude within a plausible range, 
with the constraint on {av) varying by a factor of < 10 
as B is varied from 10 - 100 /iG. (Taken alone, the radio 
and gamma rays bounds are individually quite sensitive 
to the assumed magnetic field amplitude.) Our analysis 
is distinct from previous work in this area by the use of 
a new synthesis of Galactic center radio data assembled 
in Ref. [75]. 

Our constraints on velocity-averaged DM annihilation 
cross sections are conservative as we (i) do not remove 
(known) astrophysical contributions to the radio emission 
and (a) do not include a contribution to the total radio 
emission from DM annihilation electrons and positrons 
along the line-of-sight but out of the GC. 

Despite these conservative assumptions, our results are 
at least one order of magnitude and up to two orders of 
magnitude better than the most directly comparable pre- 
vious limit (obtained by Borriello et al. for the same DM 
distributions and electron injection spectrum, but con- 
sidering all-sky radio fluxes). Moreover, our constraints 
rule out a sizable portion of the parameter space that has 
been invoked to explain the various positron anomalies. 
(Sec, e.g. Refs. [122, 124] for recent determinations of 
the PAMELA/Fermi perferred regions.) 

For an NFW profile and a xX ~^ QQ annihilation pro- 
cess, we find the allowed boost factor is < 1 at 10 GeV 
and < 100 at 1 TeV. For xx — ^ e+e~, the allowed boost 
factor is < 1 at 10 GeV and < 1000 at 1 TcV. Note that 
these constraints apply to boost factors generated via 
whatever mechanism, be it a Breit-Wigner resonance, 
the Sommerfeld effect, or an enhancement due to DM 
substructure ^. 



^ An interesting exception is the case where we can not reliably 
extrapolate between the local DM annihilation rate and that in 
the Galactic center, e.g., because of a variation in the dispersion 
velocity [125] or a local clump of dark matter. 
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Appendix A: Electron transport 

Much research concerning cosmic ray electrons adopts 
the approach of taking measured or inferred values of 
the "diffusion coefficient" (for cosmic ray electrons at a 
particular energy), assuming some energy scaling for this 
quantity, and then deriving characteristic distances over 
which electrons might diffuse over their cooling times. 
There are many uncertainties attendant upon this sta- 
tistical approach, particularly for the Galactic center en- 
vironment where evidence points to ambient magnetic 
fields stronger and more turbulent than found typically 
in the disk. We prefer to take the approach of referring to 
direct modeling (by solving the Lorentz force equation) of 
individual particle trajectories for the GC environment. 
Here we rely on the work of Wommer et al. [126] who 
consider proton propagation in this environment. From 
the work of these authors we infer that ©(TeV) energy 
protons might travel < 0.2° or 30 pc (at the GC dis- 
tance) over their radiative lifetimes for an environment 
with a 10 /LtG ambient field and uh = 100 cm~^. In this 
environment the TeV proton cooling timescalc (due to 
ionization and p — p collisions) is ~ 4 x 10^ years. In 
contrast, for GeV scale electrons (synchrotron-radiating 
at about GHz frequencies), the cooling timescalc is less 
than this at < 2 x 10'^ years for the weakest magnetic field 
(10 /uG) and most tenuous gas environment {jih — 3 
cm~^) we consider. Given that both TeV protons and 
GeV electrons are ultra-relativistic and that diffusion dis- 
tance will decline towards lower particle energies, we in- 
fer from Ref. [126] that GeV electrons will travel well 
less than 0.2° or 30 pc in their radiative lifetimes. Given 
this angular scale is well smaller than both the DNS and 
HESS fields we conclude that it is safe to conduct our 
calculations in the thick target limit. 
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